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Abstract. We present the B-band luminosity function and comoving space and luminosity densities for a sample 
of 2779 I-band selected field galaxies based on multi-color data from the CADIS survey. The sample is complete 
down to /sis = 22 without correction and with completeness correction extends to /sis = 23.0. By means of a new 
multi-color analysis the objects are classified according to their spectral energy distributions (SEDs) and their 
redshifts are determined with typical errors of 8z < 0.03. We have split our sample into four redshift bins between 
z = 0.1 and z — 1.04 and into three SED bins E-Sa,Sa-Sc and starbursting (emission line) galaxies. The evolution 
of the luminosity function is clearly differential with SED. The normalization <f>* of luminosity function for the 
E-Sa galaxies decreases towards higher redshift, and we find evidence that the comoving galaxy space density 
decreases with redshift as well. In contrast, we find <j>* and the comoving space density increasing with redshift 
for the Sa-Sc galaxies. For the starburst galaxies we find a steepening of the luminosity function at the faint end 
and their comoving space density increases with redshift. 
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1. Introduction 

The luminosity function, specifying the density of galaxies 
within a given comoving volume as a function of type and 
magnitude, is a basic ingredient in describing the galaxy 
population. A determination of the luminosity function 
at different redshifts describes in a global way the evolu- 
tion of the galaxy population with cosmic look back time. 
Luminosity function data, combined with other data such 
as number counts and direct imaging data, form the ob- 
servational basis for probing the statistical evolution of 
galaxies. 

The present epoch [z < 0.2) luminosity function is 
important as reference point of the evolution; there has 
been much controversy about its shape and normalization , 
but now it appears well established ( |Zucca et al. 1999 ) 
that the faint end Mb > — 16 of the luminosity function is 
dominated by galaxies of later morphology, later spectral 
type, bluer color and stronger line emission. 

To map galaxy evolution over the last half of the uni- 
verse's age, three surveys with samples of several hun- 
dred galaxies have been publis hed in the last years: the 
CFRS survey ([Lilly et al. 1995| ) with 591 I-band selected 
galaxies with measured redshifts extending out to z » 1, 
the autofib survey ( Ellis et al. 199C ) with > 1700 galax- 



ies with redshifts z < 0.75, and the CNOC2 survey 



(Lin et al. 1999) with over 2000 R-band selected galax- 
ies with redshifts z < 0.55. Lilly et al. 1995 found no 
evolution of the luminosity function for galaxies ear- 
lier than Hubble type Sb but a steepening of the faint- 
end luminosity function with increasing redshift for later 
galaxy types, i.e. the galaxy population is richer in dwarfs 
with increasing redshift. Such a steepening of the lumi- 



nosity function has also been found by Ellis et al. 1996. 
In contrast, the CNOC2 survey claimed to find lumi- 
nosity evolution for early type galaxy luminosity func- 
tion and density evolution for the late type galaxy lu- 
minosity function, which is not in agreement with the 



CFRS results. Kauffmann, Chariot & White, 1996 found 
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strong evolution in the population of the early type galax- 
ies by re-analyzing the CFRS data, in the sense that 
at z = 1 two-thirds of the nearby early type galaxy 
population have drop ped out of the sample. However, 
[Totani fc Yoshii 1998| did not confirm this conclusion, but 
found that the data are consistent with passive evolution. 

In this paper we present an independent determination 
of the luminosity function based on a new, large sample 
of galaxies. The CFRS, autofib and CNOC2 surveys are 
based on direct imaging data from which objects are sep- 
arated into stars and galaxies based on morphology; then 
redshifts are measured for galaxies down to a certain lim- 
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iting magnitude in a given band. As many faint galaxies 
are barely resolved, size selection effects, which can not 
be easily understood or modelled, enter in the determi- 
nation of the luminosity function. The approach we take 
here is completely different. We use multi-color data and 
determine galaxy type and redshift from comparison to a 
library of color indices. This automatic process can be sim- 
ulated and therefore we can reliably determine the com- 
pleteness of our data. A further advantage of our 12-color 
data is that we do not have to apply k-corrections by fits 
of polynomials or the like; rather we interpolate measured 
data. The disadvantage of our multi-color method are red- 
shift errors which are larger than spectroscopic ones, but 
we can show that they do not affect the luminosity func- 
tion. 

2. The Data Sample 

Since CADIS will be described in detail elsewhere 
(Meisenheimer et al., in preparation), we only briefly sum- 
marize the most pertinent information necessary for this 
paper. Once finished, CADIS will cover seven 100D' fields 
in different regions of the sky which avoid obvious fore- 
ground stars and for which the IRAS cirrus at lOO^m is 
< 2MJy/sr. Our present analysis is based on three widely 
separated fields (lh, 9h and 16h field, see table ||). 

Table 1. The centers of the fields 



field 


02000 


(52000 


1 h 


^47^3383 


2°19'55" 


9 h 


gh 13 m 46 s 5 


46°14'20" 


16 h 


16 h 24 m 32?3 


55°44'32" 



CADIS is a multi-color survey, using three broad band 
(B, R, J or K') and up to 13 medium band (S\ — 250 
to 500A) filters with central wavelengths from 3900 to 
9150A. Although these filters are tailored for detection 
of high redshift emission line objects, the resulting multi- 
color photometry allows excellent classification of the ob- 
jects according to their SEDs and determination of their 
redshifts. 

Total integration times in a narrow band filter are on 
the order of 20 ksec on the 2.2m telescope, split into several 
dithered exposures; these long exposure times result in 
10er detection limits of m > 23.0 in all filters. 

For object detection we use the SExtractor software 
package (Bertin & Arnouts 1996) on the coadded frames 
for each filter. Photometry of the detected objects is 
done with our own software package MPIAphot. This 
photometry package optimizes the signal-to-noise ratio 
by deriving the fluxes above local background with a 
weighted sum and taking seeing variations in to account 
( Roscr fc Meisenheimer 1991 , Thommcs 1996 ). Since the 
photometry is an aperture photometry with apertures 1.3 
times seeing measured on the individual frame, the mag- 
nitudes of large objects must be corrected for aperture 



losses. These corrections, which are dependent on the mor- 
phology of the galaxy, were derived from photometry on 
simulated images where the true magnitudes were known. 
Because the photometry is performed on individual rather 
than on coadded frames, realistic estimates of the photo- 
metric errors can be derived straightforwardly. Conversion 
of fluxes to magnitudes is based on observations of spec- 
trophotometric standard stars and synthetic photometry, 
i.e. integration of the instrumental response function over 
the standard star spectrum. From the loci of the mea- 
sured stellar main sequence we estimate that the system- 
atic errors in absolute photometry are less than 0.03 mag 
in all filters. The magnitudes we use are Vega magnitudes. 
For each object our photometric catalogue thus contains 
(among other data) 10 to 13 magnitude entries, their er- 
rors and 9 to 12 color indices. 

The color indices in the object catalogues are used to 
classify the objects into stars, galaxies and QSOs by com- 
paring them to a color index library derived from the SEDs 
of stars, galaxies and QSOs. Note that we do not apply 
any morphological star/galaxy separation or use other cri- 
teria, the classification is entirely spectroscopic. The stel- 
lar library is taken from Pickles (1998)[ The QSO spec- 
tral library is based on the high signal-to-noise template 
spectrum Francis et al. (1991), but also includes different 



continuum slopes and line-widths; the final QSO spectral 
library contains templates with redshifts up to z < 6. The 
spectral library for galaxies is derived from the mean aver- 



aged spectra of Kinney et al. (1996). From these, a grid of 



redshifted template spectra was formed covering redshifts 
from z = to z = 2 in steps of Sz — 0.01 and 100 dif- 
ferent intrinsic SEDs, from old populations to starbursts. 
Our final library of color indices contains entries for 131 
star-, 45150 QSO- and 20100 galaxy templates. Using the 
minimum error variance estimator, each object is assigned 
a type (star, QSO, galaxy), a redshift (if it is not classified 
as star) and an SED. The formal errors in this process de- 
pend on magnitude and type of the object and are on the 
order of 5z = 0.03 and 5SED = 2 (of 100). The dispersion 
in the estimated values does not increase with redshift for 



2000 



z < 1.04. Full details are given in Wolf et al. 

Down to the magnitude limit 771/815 > 23.0 in a 
medium band filter centered at A = 815nm, 4303 objects 
were detected, 835 of these were classified as stars. 69 as 
QSOs and 3399 as galaxies, 2779 of which are within the 
redshift range z = 0.1 and z — 1.04. For 257 objects (i.e. 
6 %) neither redshift nor SED could be estimated. 

In Fig. [l] we show the differential number counts of 
objects which are classified as galaxies and for which the 
classification algorithm can determine a redshift. As this 
figure shows, our data are effectively complete to m/815 = 
22; only in the faintest brightness bin we have to correct 
for incompleteness. The slope in the number counts from 
Tfi — 18 to m = 22 is 0.403 ± 0.0 18 and agrees well with 
the results of Nonino et al. 1999 ; fitting their data in the 
same magnitude range we obtain a slope of 0.39. Fig. || 
shows the completeness of the classification as function 
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of magnitude derived from the simulations described in 
section 4. 
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Fig. 1. Differential number counts in the 7815 filter. The 
data points include only those objects which were used in 
our analysis, i.e. objects which were classified as galaxies 
and for which a redshift was estimated. The straight line 
is the fit to these data from m/sis = 17. to m/ 8 i5 = 22.5 



Lilly et al. 1995 



have pointed out the advantages of selecting a sample 
in the I-band, which are (i) good match to B selected local 
samples since the B band is redshifted to I for z w 0.9 and 
(ii) minimization of selection effects since k-corrections are 
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Fig. 2. The completeness of the classification process de- 
rived from simulations (see sect. ||). 



much smaller and depend less on the SED of the object 
in the I-band than in shorter wavelength filter bands. Put 
another way, surveys that select the objects at shorter 
wavelengths are more affected by selection effects, since 
at given measured B-magnitude galaxies of different mor- 
phological types will contribute to the rest frame B band 
luminosity function. 

Since the redshifted B-band always falls between two 
of our filters for redshifts z < 1.04, the rest- wavelength 
B-magnitude is computed for all objects by interpolation 
between the 2 adjacent CADIS filters. We emphasize that 
we do not have to derive K-corrections from model spectra 
or polynomial extrapolations, rather we simply interpolate 
measured data. 

Figure U shows our basic observational result, the loca- 
tion of the galaxies in the Mb — z diagram. For consistency 
with earlier work the absolute B-magnitudes have been 
computed using H = 100 kms~ 1 Mpc~ 1 7 qo = 0.5. Note 
that at z — 1 the inferred absolute magnitudes would be 
0.6 mag brighter for an Hq = 100 kms~ 1 Mpc~ 1 ,£l m — 
0.3, Da = 0.7 cosmology. 

The upper right portion of this M-z diagram is devoid 
of galaxies which reflects our magnitude limit; since this 
limit is set to m = 23.0 in the Isi5 filter, it appears as a 
sharp cutoff in the Mb magnitudes only around z « 0.8. 
The region z < 0.3 is sparsely populated which is due to 
the relatively small field of view covered by our data; we 
have therefore concentrated our analysis to z > 0.3. The 
distribution of the galaxies in Fig. || shows conspicuous 
horizontal bands; all three fields combined show a more 
homogeneous distribution. Inspection of the spatial dis- 
tribution of the galaxies in our fields shows that the 9h 
field contains a sheet of galaxies near z = 0.2. Since the 
local bin contains only very few galaxies when the 9h field 
is taken out, we have therefore restricted our discussion 
to the redshift range z = 0.3 — 1. A minor contribution 
to the inhomogeneity in redshift distribution is caused by 
the special choice of filters, which favors some redshifts. 



However, this can be corrected very well (see sect. 3.2) 



3. Determination of the luminosity function 

3.1. Methodology 

The Schechtcr function is a convenient way to parametrize 
the luminosity function of galaxies. However, Fig. || illus- 
trates that the Schechter function cannot be equally well 
determined at all redshifts as the segment of the luminos- 
ity function which is covered by data changes with red- 
shift. We would expect a poor determination of M* for 
z = 0.1 — 0.3 since there are few galaxies brighter than 
M* in our relatively small field of view, and a poor de- 
termination of a for z k, 1 since there are fewer galaxies 
fainter than M* due to the magnitude limit of the survey. 
Since the parameters of the Schechter function are covari- 
ant, the parametric description of the luminosity function 
may become ill defined. Thus it is important to use both 
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Fig. 3. The Mb — z diagrams for the galaxies in our survey for the three fields individually and combined. Absolute 
magnitudes calculated for a H = 100 kms~ 1 Mpc~ 1 , qo = 0.5 cosmology. 



parametric and non-parametric estimations of the lumi- 
nosity function. 



and within a redshift interval 
the volume from 



We compute 



Several estimators to derive the luminosity function 
from a set of M-z data have been developed in the past. 
We used a modified V ma x estimator as non-parametric es- 
timate of the luminosity function and a modified maxi- 
mum likelihood estimator (STY) as parametric estimate 
(see eg. Efstathiou, Ellis & Peterson 1988). These modifi- 



cations include an explicit completeness correction, which 
is given by a function C(Zj, SEDi, M,) which is the prob- 
ability with which a galaxy at redshift 2, with SEDi and 
absolute magnitude Mi enters our sample. The determi- 



nation of C is described in Section 3.2. 



In its standard form, the V max estimator is computed 
from §{M)dM — J2i v- (m) wnere ^ i s tne total comoving 
volume in which galaxy i could lie and still be included 
in the sample. The volume is determined by the redshift 
and apparent magnitude boundaries, and the sum is taken 
over all galaxies in the [M - dM/2, M + dM/2] interval 



Vi = duj 



C{z,SED h Mi)^-dz 
ailaz 



(1) 



where dV/ (dQdz) is the comoving volume element per 
unit solid angle measured in steradians and unit redshift, 
and dto is the area in steradian covered by the survey. 
Thus the completeness function not only corrects for de- 
tection incompleteness, but also sets the limits of the in- 
tegral since C(zi, SEDi, Mi) — if an object lies outside 
of the survey criteria. The square roots of the variances 
var(Q) = J2i v{M) 2 are use d as 1°" error bars on the V max 
data points. 

In the STY method the probability pi of observing 



galaxy i is given by pidM = §{Mi)dM/ J 



Mi, 



$(M)dM 



where $(M) is the luminosity function and Mn m t z \ is the 
absolute magnitude corresponding to the magnitude limit 
of the survey for the redshift of the galaxy. The complete- 
ness correction has been incorporated by replacing <J>(Mj) 
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with $(Mi)C(zi,SEDi 
all galaxies is given by 



Mi) so the likelihood to observe 



n 

i=l 



MMjCj^SEDuMi) 
/+" ^(M)C(zi, SED U M)dM 



(2) 



where N b s is the total number of galaxies in the 
(sub)sample. Here again the completeness function C also 
sets the limits of the integral. The best parameters a 
and M* are found by maximizing InL. Errors in these 
paramaters are derived from the likelihood contours at 
InL— ^A% 2 where A% 2 is the change in \ 2 corresponding 
to the desired confidence level for 2 degrees of freedom. It 
should be noted that these errors underestimate the true 
ones, because they do not account for the cosmic inhomo- 
geneities apparent from Fig. |[ Since the normalization 
factor $* cancels out in equation ^ we have determined it 
from $* = N obs /N pred where 




Fig. 4. The completeness correction function C for galax- 
ies earlier than Sbc(left) and galaxies later than Sbc 
(right). The grey scale coding is linear from C=0 (white) 
to C=2 (black). 



obviously given by the ratio objects classified to be in that 
cell Nsed'. z'. M , t° objects initially in it: 



N, 



pred 



= duj 



+ OC 



dV 
dQdz 



$(M)C(z,M)dzdM (3) C(z l ,SED l ,M l ) = 



NsED'. 



NsEDi,Zi,M E 



(4) 



is the predicted number of galaxies using the derived 
parameters a and M* in the luminosity function $(M), 
and C(z, M) is the mean of C over all SEDs. The error of 
$* was determined through bootstrapping. 

3.2. Completeness Correction 

The completeness function C was determined as follows. 
For a given set of colors and their errors, the color 
classification algorithm (Wolf et al. 2000) estimates red- 
shift and SED for each object. After this classifying pro- 
cess an object drawn from a given point (SEDi, 2j, MgJ 
in the 3 dimensional parameter space will appear at 
(SEDi, z\, M' B ). Even for random measurement errors the 
expectation value of this latter point in the parameter 
space is not necessarily identical to the original one, so 
the observed parameter space may have systematic over- 
or underdensities introduced by the classification process. 
Furthermore, the object may not be classified at all which 
indeed sets the completeness limit to our data. However, 
both completeness limit and possible classification system- 
atics can be effectively calculated and corrected by Monte 
Carlo simulations of the observations. 

From the galaxy template library we have constructed 
a grid of 625191 input "objects" brighter than rn/815 = 
23.5 in the 3 dimensional parameter space SED — z — Mb 
where < SED < 100, < z < 1.1 , -12 < M B < -24. 
Apparent magnitudes were calculated for each object cor- 
responding to its redshift and Gaussian noise was added, 
matching the real measurement errors. Then this cata- 
logue of simulated objects was passed through the classify- 
ing process and a new Mb', using the redshift z[ estimated 
by the classifying routine was calculated. 

Dividing the SED — z — Mb parameter space into cells, 
the completeness function C(z, SED, Mb) in each cell is 



The parameter space was devided into 25500 cells. 

The determination of C was done separately for the 
two cosmologies; in Fig. [| we show 2 D averages for 
SED = 0-50 and SED = 50 - 100 for the Q m = 
1,0a = cosmology. It is obvious that the classification 
process favors some regions of the redshift space and dis- 
criminates against others; this is probably caused by the 
special choice of the wavelengths of the filters we are using. 
However, since our definition of C allows C > 1, this ef- 
fect is fully taken into account in the corrections for V max 
(Equation pi) and in the maximum likelihood estimation 
(Equation g) . 

4. Monte Carlo Simulations 

The effect of photometric errors on the luminosity func- 
tions was studied by Monte Carlo simulations. An advan- 
tage of our purely photometric classification is the possi- 
bility to simulate the complete process, including the red- 
shift assignments to the objects. Though the errors of our 
photometric redshifts are on the order of Sz = 0.04 and 
thus larger than typical spectroscopic redshift errors, one 
would expect that their effect is small since they lead to 
errors in absolute magnitude SM = 0.08 which are much 
less than the 1-mag bins we are using for calculating the 
luminosity function. Nevertheless, we have examined the 
effect of redshift errors by numerical simulations. 

Simulated galaxies were distributed in the M-z plane 
with a redshift distribution that corresponds to a constant 
comoving volume density and with absolute magnitudes 
drawn from a Schechter function with M* = —19.5, a 

1 .0 for all types of galaxies. The SEDs were randomly 
distributed in the interval [0,100]. Apparent magnitudes 
were calculated using the fl m = 1, Oa = cosmology and 
noise added to the color indices according to the actual 
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measurement process for all artificial galaxies. Th 
jects were then run through the classifying proce 
to/815 > 22.0 94% of the input galaxies were class 
such, and 86% of the galaxies have redshift errors < 
In the center of our last magnitude bin m/gi5 : 
the corresponding numbers are 82% and 72%, respe 
The classification completeness is shown in Fig. |2|. 
manner we have formed 5 samples of galaxies w 
number of galaxies in each sample roughly as in ox. 
For these samples we have then determined the lurr 
function as for the real galaxies. 

The agreement between the input and the re< 
Schechter parameters is excellent in the redshift 
z = 0.3 — 0.9 where we find displacements in 
teristic magnitude \SM*\ < 0.05 and slope \8a\ ; 
As expected, the parameters are less well recove 
z = 0.1 - 0.3 and z > 0.9: we obtained M* = 
a = -0.82 for z = 0.1-0.3 and M* = -19.84, a = 
for z > 0.9. The numbers given here apply to bot 
and late type galaxies; we found no significant depe 
on SED type. 

These simulations thus show that we can reco 
luminosity function parameters very well for z = 0. 
The displacements of the recovered parameters are 
than their errors in the actual data. For redshifts 
the z = 0.3 — 0.9 range the description of the k 
ity function by the parameters of the Schechter fi 
becomes poor. However, even if the parameters a 
determined, it may be dangerous to draw conclus 
the luminosity function from the parameters alone 1 

they are coupled: for example, a decrease in M* L 

pensated by an increase in a - the error ellipses in Figs. 
^ to [ll] demonstrate this very clearly - and consequently 
by a decrease in <jf . Therefore the non-parametric form of 
the luminosity function is indispensable. 

5. Results 

5.1. The local luminosity function 

Because of the sheet of galaxies contained in the 9h field 
and the resulting problems (see sect. our results for 
the local redshift bin are preliminary until CADIS is 
completed. Nevertheless we show the luminosity function 
for all galaxies with z = 0.1 — 0.3 in Fig. |j. The de- 
rived Schechter function is described by <f>* = 0.0156 ± 
0.004, M* = -19.2 ±0.27, a = -1.30 ± 0.069. This is 
agrees within the errors with an eyefittcd mean of re- 



cent determinations by 


Zucca et al. 1999 




Loveday 1997, 


Marzkc & da Costa 1997 




Marzkc, Huchra & Gellcr 1994 


and 


Folkcs et al. 1999 


which is described by (f>* — 



0.017, M* = -19.6, a = -1.15. 



5.2. Evolution of the luminosity function 

It has been well established that the local luminosity func- 
tion is different for early and late type galaxies in the sense 
that the slope a « —0.8 for early type galaxies which 
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Fig. 5. The luminosity function of the 439 galaxies in our 
sample with z = 0.1 — 0.3. The data points were derived 
from the V max method, the solid line is the Schechter func- 
tion derived from the STY method, the dashed line is 
the mean of recent determinations of the local luminosity 
function. 



are dominated by old populations and a « —1.8 for late 
type galaxies which are dominated by starforming pop- 
ulations (Marzke et al. 1994), i.e. the starforming galaxy 
population is richer in dwarfs. Furthermore, other redshift 
surveys have found a dependence of the evolution of the 
luminosity function on galaxy type. We have therefore di- 
vided our sample into three SED bins, E-Sa galaxies with 
SED < 30, Sa-Sc galaxies with 30 < SEP < 75 and 
starforming galaxies with SED > 75 (|Wolf et al. 2000, 



Kinney et al. (1996) ). We have further divided our sam- 
ple into three redshift bins z — 0.3 — 0.5, z — 0.5 — 0.75 
and z = 0.75 - 1.04. 

To facilitate comparison with other work we 
have derived the luminosity function for an Hq — 
100 kms~ 1 Mpc~ 1 , qo — 0.5 cosmology. We also give the 
results for a cosmology with TO = 0.3, f^A = 0.7; in order 
to show the effects of cosmic acceleration in a transpar- 
ent manner, we have again used Ho = 100 kms Mpc . 
For Hq = ^2kms~ 1 M'pc~ 1 the conversion relations are 



M, 



62 



Mi 



-1.02 and, 



'62 



^100 



x4.2. Since the effects of 



the cosmology on the luminosity functions and comoving 
space densities are mild, we concentrate the discussions 
on the go = 0.5 cosmology. The results are summarized in 
Figs. | to m and table |72. 

As Figs. U - ^ show, the luminosity functions of 
early type galaxies E-Sa have indistinguishable M* and a 
whereas <fi* decreases with redshift. For the Sa-Sc galaxies 
M* and a are constant, too, but <fi* is increasing with red- 
shift. Note that the luminosity function of the starburst- 
ing galaxies with z = 0.75 — 1.04 depends sensitively on 
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a few bright galaxies in the M = —22 bin. Nevertheless, 
a steepening of the luminosity function with redshift is 
clearly evident. 

To avoid the covariances in the parametric form of the 
luminosity function, we have also analysed the comoving 
galaxy space density and luminosity density as function 
of redshift in a way which is completely independent of 
the luminosity function. Since our data are virtually com- 
plete to Mb = —18.5 (see Fig.||) for all redshifts, we have 
integrated the numbers and luminosities, respectively, for 
galaxies brighter than this limit without completeness cor- 
rection. The resulting comoving densities are shown as 
function of redshift in the Figs. |l2] and |l3] and the param- 
eters of linear fits to these data are given in table |5.2| . 

For the qo = 0.5 cosmology, we find a decreasing 
comoving space density of the early type galaxies. The 
space density increases with redshift for the Sa-Sc galax- 
ies as well for the starburst galaxies. The behaviour of 
the B-band comoving luminosity density is similar. For an 
fl m = 0.3, ft\ = 0.3 cosmology, the trends in the data are 
similar but with reduced gradients. 

As customary, we have calculated the errors based on 
Poisson statistics. This is correct as long as the total vol- 
ume sampled by the survey is considered, where large- 
scale structure averages out. The scatter in figjl] agrees 
with Poisson statistic. However, if the survey is devided 
into smaller subvolumes, cosmic variance due to cluster- 
ing of galaxies increases the actual uncertainties. Fig.[l4| 
shows the field to field variations using the same SED and 
redshift bins as applied in the determination of the lumi- 
nosity function. The scatter in figsjlj and |l3| is higher than 
Poisson by factors of a few. If this is taken into account, 
the decrease in comoving galaxy space density for the E- 
Sa galaxies with increasing redshift is statistically not sig- 
nificant, the increasing space density of the Sa-star burst 
galaxies towards earlier epochs, however, still is signifi- 
cant. We emphasize that this problem of clustering exists 
with all previous surveys; quite contrary, our survey is 
much larger and will be still larger once all targeted fields 
are available. 



5.3. Comparison with other work 

Since the luminosity function depends strongly on galaxy 
type, a comparison between results obtained from sam- 
ples selected with different selection criteria is difficult if 
not impossible. In a similar way, division into differing 
subsamples causes differences in the luminosity function. 
Furthermore, different ways to estimate K-corrections can 
introduce redshift dependent effects. Photometric errors 
should be negligible if CCD data are used. 



The CFRS ( [Lilly et al. 1995[ ) is the survey which com- 
pares most directly to CADIS: their objects are also I- 
band selected across an identical redshift range. Since the 
number of galaxies included in our analysis is nearly four 
times as large, we have better statistics. Furthermore, our 
survey goes deeper by 0.5mag and so the range in lu- 



minosities covered is larger, which enables better fits of 
the Schechter function. When we divide our sample in the 
same two SED bins as CFRS, our findings agree very well 
with their main conclusions, namely no evolution of the 
early type galaxies (E-Sbc) and brightening and steepen- 
ing of the luminosity function with redshift for later types. 
Inspection of the data points shows that the two surveys 
agree within the errors (though our errors are smaller). 



Lilly et al. 1996 



have determined the evolution of the luminosity den- 
sity from the CFRS data. For the z — 0.35 bin they ob- 
tained logL = 19.48Vt / iJz _1 (the value we use here is 
their 'directly observed' 4400 value, converted to Hq = 
100). Using Lq(B) = 3.4 lO 11 ^^" 1 we obtain log L = 
19.38 WHz' 1 at z = 0.35 from the fits in table |J 
value includes only galaxies brighter than Mb = - 



This 
18.5; 



the correction factor to the full population down to Mb = 
— 14 determined in the two lowest redshift bins is a factor 
of 2, resulting in logL = \9.Q8WHz~ 1 which agrees wit hin 
the errors with the value obtained by Lilly et al. 1996. At 



z — 1 Lilly et al. 1996 find an increase in comoving lumi- 
nosity density of the blue galaxies by a factor of 2.9 over 
z — 0.35 compared to our data which give a factor of 3.4 
for the late types, again well within the errors. We there- 
fore conclude that our deeper and larger redshift survey 
is in excellent agreement with the CFRS. 



The CNOC2 Lin et al. 1999 sample is R-band selected 



and complete to R = 21.5. It contains about 2000 galaxies 
with redshifts 0.12z < 0.55. Their luminosity functions do 



not agree well with those of the CFRS (see Lin et al. 1999) 
and hence also not with ours. 

Possible explanations for differences to CFRS have 
been discussed by Lin et al. 1999 but the cause remains 
unclear. It seems plausible that sample selection effects 
and k-corrections do at least partly cause the differences, 
since the CNOC2 sample is R-band selected whereas 
CFRS and CADIS are I-band selected. 

The interpretation of the CNOC2 data is based on 
the evolution of the parameters of the Schechter function; 
since these are coupled, the parameters describing the evo- 
lution are coupled, too, which may give misleading results 
(see below). Also, Lin et al. 1999 keep a fixed, although 



the luminosity functions of the late type (starburst) galax- 
ies clearly show steepening. One should further note that 
the CNOC2 results are obtained over a smaller redshift 
range than CFRS or CADIS. Additionally, the errors of 
the luminosity function determinations in all available sur- 
veys are large and the significance of the results is low. 



The autofib survey ( [Ellis et al. 1996[ ) is a blue selected 
sample and contains redshifts for « 1700 galaxies. A steep- 
ening of the luminosity function from redshifts z = 0.02 
to z — 0.75 has been claimed, although the error el- 
lip ses for the lowe st and highest redshifts bins (Fig. 11 
in Ellis et al. 1996| ) overlap on the 3<r level. 
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Table 2. Parameters of Schcchtcr function fits to the data. 







— 1 












sample 


z range 










M* 


a 


E-Sa 


0.3-0.5 


0.901 10" 


•2 


± 0.559 10" 




-19.287 ± 0.219 


0.187 ± 0.224 


E-Sa 


0.5-0.75 


0.746 1(T 


2 


±0.48710" 


3 


— 19.355 ± 0.180 


—0.25 10 ± 0.194 


E-Sa 


0.75-1.04 


0.455 1(T 


2 


±0.562 10" 


3 


-19.360 ± 0.337 


0.282 ± 0.52 


Sa-Sc 


0.3-0.5 


0.13710" 


1 


± 0.355 10" 


2 


-19.660 ± 0.295 


-0.815 ± 0.129 


Sa-Sc 


0.5-0.75 


0.159 10" 


1 


±0.216 10" 


2 


-19.335 ± 0.158 


-0.615 ± 0.128 


Sa-Sc 


0.75-1.04 


0.182 10" 


1 


± 0.462 10" 


2 


-19.770 ± 0.268 


-0.807 ± 0.295 


starburst 


0.3-0.5 


0.17710" 


1 


±0.71510" 


2 


—18.690 ± 0.292 


-1.282 ± 0.209 


starburst 


0.5-0.75 


0.281 10" 


2 


± 0.273 10" 


2 


—20.27 ± 0.400 


-1.857 ± 0.168 


starburst 


0.75-1.04 


0.265 10" 


1 


±0.102 10" 


I 


-19.441 ± 0.227 


-1.727 ± 0.271 






Q m = 0.3, ST2 A = 0.7 








sample 


z range 


<f>* 








M* 


a 


E-Sa 


0.3-0.5 


0.488 10" 


•i 


±0.275 10" 


3 


-19.602 ±0.178 


0.207 ±0.178 


E-Sa 


0.5-0.75 


0.31510" 


1 


± 0.332 10" 


3 


-19.871 ±0.276 


-0.475 10 _1 ±0.221 


E-Sa 


0.75-1.04 


0.144 10" 


2 


±0.30110" 


3 


-19.705 ±0.316 


0.625 ± 0.587 


Sa-Sc 


0.3-0.5 


0.831 10" 


2 


±0.139 10" 


2 


-19.911 ±0.193 


-0.782 ± 0.892 10" 1 


Sa-Sc 


0.5-0.75 


0.69710" 


2 


±0.971 10" 


3 


-19.832 ±0.159 


-0.605 ±0.159 


Sa-Sc 


0.75-1.04 


0.676 10" 


2 


±0.11510" 


2 


-20.232 ±0.194 


-0.667 ±0.219 


starburst 


0.3-0.5 


0.11910" 


1 


±0.360 10" 


2 


-18.970 ±0.233 


-1.141 ±0.184 


starburst 


0.5-0.75 


0.211 10" 


2 


±0.162 10" 


2 


-20.365 ± 0.299 


-1.730 ±0.147 


starburst 


0.75-1.04 


0.108 10" 


1 


± 0.329 10" 


2 


-19.892 ±0.203 


-1.480 ±0.276 



Table 3. Parameters of linear fits of the form a ± b(l ± z) to the comoving galaxy density and B-band luminosity 
density. Absolute terms a are given in the 2., slopes b in the 3. column 



fi m = l,n A = 0,Q A = 

a b 



number density E-Sa 


0.669 10"* 


± 0.192 10 -2 


-0.174 1Q- 2 ±0.111 10" 


•i 


number density Sa-Sc 


0.681 10" 3 


±0.339 10" 2 


0.744 10" 


2 ±0.205 10 -2 




number density starburst 


-0.21610" 


1 ±0.191 10" 2 


0.184 10" 


1 ±0.129 10" 2 




luminosity density E-Sa 


-0.138 10 s 


±0.163 10 8 


0.354 10 s 


±0.103 10 s 




luminosity density Sa-Sc 


-0.153 10 9 


± 0.300 10 8 


0.179 10 9 


±0.193 10 s 




luminosity density starburst 


-0.198 10 9 


±0.163 10 s 


0.167 10 9 


±0.113 10 8 






f2m = 0.3, 
a 


ft A = 0.7 


b 






number density E-Sa 


0.630 10 -2 


±0.104 W- 1 


-0.28910-^0.583 10" 


3 


number density Sa-Sc 


0.871 10 -2 


±0.181 10~ 2 


-0.175 10~ 2 ±0.105 10" 


2 


number density starburst 


0.630 10" 2 


±0.104 10" 2 


-0.269 10~ 2 ±0.583 10" 


3 


luminosity density E-Sa 


0.309 10 8 ± 0.100 10 8 


-0.640 10 7 ± 0.582 10 7 




luminosity density Sa-Sc 


0.899 10 7 ± 0.164 10 7 


0.31710 s 


± 0.989 10 7 




luminosity density starburst 


-0.756 10 s 


± 0.878 10 7 


0.659 10 s 


± 0.567 10 7 





6. Discussion 

We have measured the evolution of the luminosity func- 
tion; its traditional interpretation in terms of density vs. 
luminosity evolution is problematic in the context of a 
hierarchical picture. For the early type galaxies E-Sa, 
both (j>* and the comoving space density appear to de- 
crease with redshift, albeit statistically significant only 
if clustering of galaxies is neglected. Since M* and a 
are constant, this result implies that there are fewer 
early type galaxies at higher redshift. This is exactly 
what one would expect in hierarchical clustering scenar- 
ios which require that massive objects form late by merg- 
ing of smaller objects. Opposed to this scenario is the 



early formation of galaxies at high redshift in a collapse 
and a single burst of star formation followed by passive 
evolution. Observational evidence has been found for ei- 
ther scenario (see [Bchade et al. 1999| for a detailed dis- 
cussion). From our data, the E-Sa galaxy space density 
at z = 1 is lower by a factor of 1.6 (qo = 0.5) com- 
pared to z = 0. This is close to the factor of 2-3 pre- 



dicted by hierarchical cluster ing models (Kauffmann 1996, 
Baugh, Cole k Frcnck 199(1 . The increase of the space 



density of the Sa-Sc and starburst galaxies with redshift 
fits into this picture as well. 

Though our data clearly show evolution of the galaxy 
population, the interpretation of the data is not straight- 
forward. The luminosity function describes the galaxy 
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population in a statistical way and contains no direct in- 
formation on the evolution of individual galaxies. A galaxy 
may disappear from a given class (or SED bin) either be- 
cause it disappears (e.g. through merging) or because it 
changes colors (e.g. by starburst) and therefore appears in 
another bin. Therefore data like the ones presented here 
alone do not allow to distinguish between luminosity and 
density evolution. Furthermore, setting a magnitude limit 
in one pass-band may introduce systematic redshift de- 
pendent effects. Hence, modelling of different evolutionary 
scenarios including the detection and classification process 
of the survey is required. In a forthcoming paper we will 
compare our results with hierarchical scmianalytical mod- 
els of the formation and evolution of galaxies. 

For all disk galaxies (Sa-Sc and starburst) our 
data show clear evolution. It is already well estab- 
lished that disk galaxies evolve to the present day. 
Late type disk galaxies have present-day star forma- 
tion rates comparable to those at earlier cosmic times, 
early type disks have formed most sta rs at earlier times 
( Kennicutt, Tamblyn fc Congdon 1994 ). Brightening of 
gal actic disks by w 2mag at z = 1 has been found 
by [Bchade et al. 1996 which they i nterpret as luminosity 



evolution. Mao, Mo fc White 1998 have additionally used 



kinematic data and interpreted the disk brightening as de- 
crease in disk scale length ~ (1 + z) -1 . We note that the 
purely photometric classification of objects used in CADIS 
is less susceptible to systematic redshift dependent effects 
than surveys which use morphological criteria. However, 
only detailed modelling of the detection process can show 
how selection effects affect the data. 

HST imaging data are now available for a fair num- 
ber of galaxies drawn from the C FRS and LDSS 
samples with know n redshifts ( Lilly et al. 1998 , 
Brinchmann et al. 1999, Bchade et al. 1999, 



Le Fevre et al. 2000). A direct comparison with these 



studies is difficult because there is only a loose correlation 
between our spectroscopic and their morphological classi- 
fication. However, almost certainly a major fraction of our 
starbursting galaxies would be morphologically classified 
as irregular and / or 
have found an 



merging, 
increase in the 



Brinchmann et al. 1999 



proportion of irregular 



galaxies and Le Fevre et al. 2000 found an increase in 



the merger fraction with redshift. So the evolution of the 
starbursting galaxy population with redshift is supported 
by the imaging results. 



7. Conclusions 

We have determined the luminosity function of 2779 field 
galaxies based on multi-color data from the CADIS survey. 
The data base contains 9-12 color indices from optical to 
IR wavelengths for all objects from which a spectroscopic 
classification according to the SED and a redshift are de- 
rived by comparison with libraries for stellar, galaxy and 
QSO spectra. Our data are complete to igis — 23.0 and 
extend to redshifts z = 1.04. Our main findings are: 



1. The evolution of the B-band luminosity function is 
clearly differential. We find the normalization <j>* of the 
early type E-Sa galaxy luminosity function and the inte- 
grated comoving space density to be decreasing with in- 
creasing redshift although the latter with less significance. 
The normalization for the Sa-Sc galaxy luminosity func- 
tion increases with redshift as well as the comoving space 
density. The luminosity function of the starburst galaxies 
steepens towards the faint end and their comoving space 
density increases with redshift. 

2. The survey most directly comparable with our data 
base is the CFRS. When we divide or sample in the same 
galaxy types, the results agree very well. The two other 
samples which have studied the evolution of the luminos- 
ity function are not directly comparable to our data; we 
confirm the steepening of the luminosity function found in 
the autofib survey, but not all evolutionary effects claimed 
by the CNOC2 survey. 

3. The density evolution of the early and late type 
galaxy population apparent in our data is suggestive of 
merging. However, since merging increases the star for- 
mation rate and consequently the B-band luminosity, the 
interpretation is not straightforward and requires compar- 
ison to models. 
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Fig. 6. Left panel: The luminosity function of early type galaxies E-Sa for a Ho = 100 kms~ 1 Mpc~ 1 , qo = 0.5 
cosmology. V max data points with la error bars and Schechter functions derived from the STY method are shown. 
The numbers are the numbers of galaxies in the corresponding redshift bin. Right panel: la and 2a error contours. 




Fig. 7. Same as Fig. || but for galaxies of types Sa-Sc. 
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Fig. 8. Same as Fig. [j] but for starbursting galaxies 




Fig. 9. Same as Fig. || but for a Hq = 100 kms 1 Mpc 1 , fi m = 0.3, fl\ = 0.7 cosmology. 
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Fig. 11. Same as Fig. || but for starbursting galaxies 
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Fig. 12. The evolution of the comoving number density of galaxies Mb < —18.5 for the two cosmologies. Note the 
different scaling of the ordinates. 
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Fig. 13. The evolution of the comoving B-band luminosity density of galaxies Mb < —18.5 for the two cosmologies. 
Note the different scaling of the ordinates. 
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Fig. 14. The field to field variations of the number of galaxies used in the determination of the luminosity functions. 



